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We construct both normal and anomalous deterministic biased diffusions to obtain the Einstein 
relation for their time-averaged transport coefficients. We find that the difference of the generalized 
Lyapunov exponent between biased and unbiased deterministic diffusions is related to the ensemble- 
averaged velocity. By Hopf 's ergodic theorem, the ratios between the time-averaged velocity and the 
Lyapunov exponent for single trajectories converge to a universal constant, which is proportional to 
the strength of the bias. We confirm this theory using numerical simulations. 



Introduction. — Intrinsic randomness of time-averaged 
observables is of great interest in non-equilibrium statis- 
tical mechanics. Equilibrium systems are characterized 
by a few macroscopic observables, which are the time 
averages of microscopic observation functions. Macro- 
scopic observables fluctuate around their ensemble aver- 
ages in an equilibrium state. This means that the time 
averages of microscopic observation functions converge 
to their ensemble averages. In contrast, macroscopic ob- 
servables such as the diffusion coefficient and fluorescence 
intensity cannot converge to a constant but show large 
fluctuations in non-equilibrium (non-steady) phenomena 
such as anomalous diffusions and intermittent phenom- 
ena In particular, diffusion coefficients of biological 
molecules in cells [jl Q and the fluorescence times in sin- 
gle nanocrystals [3| show large fluctuations, indicating 
ergodicity breaking. 

Such intrinsic randomness of timc-avcraged observ- 
ables is known to be universal in models characterized 
by power-law trapping-time distributions with divergence 
means ^ typical example is a continuous time 

random walk (CTRW), which is a random walk with a 
random continuous trapping-time. In fact, CTRWs with 
infinite mean trapping-times show intrinsic randomness 
of the diffusion coefficient @, 0. Although ergodicity, 
i.e., time average being equal to the ensemble average, 
does not hold in such systems, time averages converge 
in distribution. This phenomenon is called distributional 
ergodicity flC^D 

In dynamical systems, distributional ergodicity is 
known to be the ergodicity in infinite measure systems 
because dynamical systems relating to stochastic models 
with infinite mean trapping-times have infinite invariant 
measures fllj ]. In infinite measure systems, the time aver- 
age for an observation function converges in distribution. 
The distribution of the time average is determined by 
properties of the observation function In particular, 
the distribution of the time average of an i 1 (/x) function 
/(x) converges to the Mittag-LefHer distribution 12j: 



n— 1 



(1) 



able with a Mittag-LefHer distribution of order a, and =>• 
denotes the convergence in distribution. The sequence 
a„ is called the return sequence, which is relevant to 
non-stationarity. In deterministic subdiffusion, where the 
mean square displacement grows sublinearly (x^) cx n Q , 
the scaling of a n is the same as a n [|| . 

In diffusion processes, an external force F generates 
a drift V. In general biased random walks including 
anomalous diffusion, the Einstein relation, 



V = 



FD 
2k B T 



(2) 



fc=0 



provided that J fd[i ^ 0, where M a is a random vari- 



holds [13|], where D is the diffusion coefficient under no 
bias, ks is the Boltzmann constant, and T is the temper- 
ature. Diffusion properties in hyperbolic chaotic dynam- 
ical systems are well known. Using the escape rate for- 
mula, one can show that the diffusion coefficient is equal 
to the difference between the Lyapunov exponent and 
Kolmogorov- Sinai entropy [l4j . Moreover, in a Lorentz 
gas, the largest Lyapunov exponent is suppressed by an 
external field and there is a relation between the Lya- 
punov exponents and a bias [IU Il6j . However, little is 
known about the relation between dynamical instability 
and the transport property under a bias . 

In this Letter, we derive a dynamical system for ar- 
bitrary biased and unbiased CTRWs. Furthermore, we 
explain the linear response to a small bias and the Ein- 
stein relation in general deterministic diffusions. Using 
the Lyapunov exponent, we obtain a relation between 
the macroscopic transport coefficient and the microscopic 
chaos in a deterministic subdiffusion. In particular, we 
relate the ensemble-averaged velocity to the difference of 
dynamical instabilities between an unbiased and a biased 
dynamical system. Using Hopf's ratio ergodic theorem, 
we show that the ratio between the time-averaged veloc- 
ity and the Lyapunov exponent converges to a universal 
constant. Moreover, we find that the universal constant 
is proportional to a bias, and the proportional constant is 
determined by the diffusion coefficient and the Lyapunov 
exponent under no bias. 

Deterministic biased diffusion model. — Dynamical sys- 
tems exhibiting diffusion are represented by one- 
dimensional map T(x) with translational symmetry, 



2 



0.5 



-0.5 



a 1 d 

Hi F 






9 


/ ti 


A' 













-0.5 -0.25 0.25 0.5 
X 



FIG. 1. (Color online) Reduced map R(x) and the corre- 
sponding original map T(x) [z — 3.0 and c = 0.3 in Eq. I|17[)]. 



T(x + L) =T(x) + L (L = 0,±1,...) 0. ThemapT(i) 
can be reduced to a map on unit interval [—1,2,1/2], 
R(x). Here we construct a dynamical system correspond- 
ing to the CTRW. The CTRW is defined by a jump 
length distribution l(x) and a waiting time distribution 
For simplicity, a jump length distribution is set to 
l(x) = pS(x — 1) + qS(x + 1) (p + q = 1). Since wait- 
ing times in CTRWs are independently and identically 
distributed random variables, we can construct a one- 
dimcnsional piecewise linear intermittent map R(x) with 
the same waiting time distribution -0(t) [l~8[. The prob- 
ability for the right jump p is represented by the length 
of the interval [c, 1/2], which is mapped to the right side 
in a neighboring cell. As shown in Fig. 1, the reduced 
map on [1/4,1/2] is given by an asymmetric piecewise 
linear map. We note that the scaling exponents of the 
waiting time distributions for the right and left jumps 
are the same. However, the waiting time distribution for 
the interval [0, 1/4] are slightly different from that for the 
interval [—1/4,0], since the derivatives R(x) at x = c + 
and c — are different in biased models [19( . When the 
mean waiting time diverges, tp(t) oc t~ l ~ a (0 < a < 1), 
R(x) has an infinite invariant measure [111 and the map 
T(x) causes subdiffusion, i.e., (x 2 ) oc n a [8|. 

Einstein relation. — The Einstein relation refers to the 
relation between the velocity under a bias and the diffu- 
sion coefficient under no bias. We obtain the Einstein re- 
lation for the time-averaged transport coefficients in nor- 
mal and anomalous deterministic diffusions. The time- 
averaged mean square displacement (TAMSD) and the 
time-averaged drift (TAD) arc defined by 



N-l 



Sx^N; x ) = — i x k+ m - xk) 



and 



fe=0 



N-l 



5x m (N; x ) = — (xic+m - Xk), 



(3) 



(4) 



k=Q 



displacement, TAMSD shows normal diffusion even when 
a < 1 Therefore, the time-averaged diffusion coeffi- 
cient is defined by D = 5x 2 n {N; x^j/m. In finite measure 
cases, the time-averaged diffusion coefficient is equal to 
the ensemble-averaged one. The ensemble-averaged drift 
(6x m ) F is given by (Sx rn ) F = (x m - x Q ) F , where (-) F 
refers to the ensemble average under a bias (c ^ 0.375). 
Let N m be the total number of jumps until time m, then 
we have (Sx m ) F = (p — q)(N m ) F . Here we assume that 
the injection to the set [-1/2, -1/4) U (1/4, 1/2] is uni- 
form (Assumption A) . This assumption is exact when the 
map is a piecewise linear map. It follows that p = 2 — 4c 
and q = 4c — 1. In normal diffusion, (N m ) F is given by 
fi(J c )m, where fi is an invariant measure of R(x) and 
J = [-1/4,1/4]. In subdiffusion, (N m ) F oc m a . Then, 
we have 



N-l 



{5x m (N)) F = (p~q)J2 



(N k+m ) - (N k ) „ - 

— = \V) F m, 



k=0 



N 



(5) 



W e- F / 2kBT ,q = W e +F / 2kBT , and An(fc) be 



where we define the time-averaged velocity by V 
5x m (N)/m 
Let p 

the total number of jumps in the interval [k, k + to], so 
that (xk+m — Xk) F — (j> — q)(An(k)) F . Here we assume 
(An(k)) F = K(An(k)) (Assumption B). K does not de- 
pend on c for a small bias because the scaling of (An(k)) F 
is the same for all c. It follows that {xk+ m — Xk) F = 
K{p - q)(An{k)) a* K(p - q)((x k+m - x k ) 2 ). We have 



where s 



(V) F = K(D)e, (6) 
8c. Thus, we obtain the Einstein relation, 
(D)F' 



(V)i 



2k B T 



F^0, 



(7) 



where F' = KF. which is not equal to F. It seems that 
the Einstein relation is violated. However, the change 
in the waiting time distribution due to the bias results 
in the violation of the Einstein relation. Therefore, the 
Einstein relation holds under the bias F', because the 
difference between the waiting time distributions for the 
right and the left does not generate a drift. We note 
that the ensemble average in the Einstein relation is 
essential in anomalous subdiffusion because the time- 
averaged transport coefficients D and V are intrinsically 
random. When the time-averaged velocity V and the 
time-averaged diffusion coefficient D are measured inde- 
pendently, V/D becomes random, namely, the effective 
temperature, FDj^ksV , becomes random j20j |. 

Relation between the Lyapunov exponent and the ve- 
locity for single trajectories. — The Lyapunov exponent is 
defined by 



N-l 



respectively. Unlike the ensemble-averaged mean square 



A(A^ )^E ln l T '^)l 



(8) 



fc=0 
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By Hopf 's crgodic theorem [Tj 



J2k=0 fm(T k x) 



Jo /2 f m dn 



^ Cm — 



Io /2 



(9) 



holds for almost all initial points of x, where f m (x) = 
T m (x) — x and g(x) = ln|T'(x)|. Therefore, the time- 
averaged velocity V(N;x) = Sx m (N,x)/m and the Lya- 
punov exponent A(iV; x) satisfy the following relation, 
V(N; x)/X(N; x) — > C m /m as N — > oo for almost all x. 
Since Sx m (N) oc to for a large m, C rn /m does not depend 
on to. Let C m /m equal xv, then we can obtain the re- 
lation between the velocity and the Lyapunov exponent 
for almost all initial points of x: 



V(N;x) 
X(N;x) 



— > xv as N — > oo. 



(10) 



In general, the ensemble average of the Lyapunov ex- 
ponent decreases with increase in the ensemble-averaged 
velocity (see Fig. 3). Nevertheless, the time-averaged ve- 
locity is proportional to the Lyapunov exponent. The 
constant xv is a universal constant, which does not de- 
pend on an initial point. Using the ensemble average and 
Eq. ([6]), we obtain 



Xv 



(V) F _ K(D) : 

(\) F (\) F 1 



fll) 



Difference of the Lyapunov exponent. — We consider 
the difference of the generalized Lyapunov exponent be- 
tween an unbiased and a biased dynamical system. The 
generalized Lyapunov exponent is the ensemble average 
of the normalized Lyapunov exponent 111 |21| 



A = ( lim -Vjfii) 



(12) 



where (•) represents the ensemble average of the ini- 
tial points. Note that the density is absolutely contin- 
uous with respect to the Lebesgue measure and satis- 
fies the condition (g(xo)) < oo. The assumption B, 
(An(k))F = K(An(k)), means that the statistical prop- 
erty of the reinjection to the indifferent fixed points is 
almost the same. That is, the generalized Lyapunov ex- 
ponent restricted to J = [—1/4, 1/4] is almost the same. 
The difference of the generalized Lyapunov exponent, 
AA, is defined by AA = <A»> - (A») F . 
By assumption A and B, 

_j n — 1 

AA = — V;«ff(a*)lj„(a*)> - (g c (x k )l Jc (x)} F ). (13) 



fe=0 



where AAjc is the difference of the generalized Lya- 
punov exponent restricted to J c . Since the injection 
to J c is uniform (Assumption A), we have AAjc = 
ln2 + phip + (1 — p) ln(l — p). The generalized velocity 
and its maximum are defined as V = lim fl _ >00 ((Jx n ) f /n a 
and Vmax = lim Il ^ 00 (A Il ) F /n a , respectively. We use 
= p — q. Then, the probability p is written by V 
m ax : p = (1 + VyVmax) /2. Using V and V max , we 



v 

and V T 



obtain the relation between the difference of the general- 
ized Lyapunov exponent and the generalized velocity: 



AA 



S 



V 



(15) 



where S(x) = (1+x) ln(l+x)+(l— x) ln(l— x). Therefore, 
the constant xv is given by 



(16) 



Note that the proportional constant is determined by the 
diffusion coefficient and the Lyapunov exponent under no 
bias. 

Example for deterministic subdiffusion. — We demon- 
strate numerical results for deterministic subdiffusion. A 
piecewise linear map is a good approximation for an in- 
termittent map. We consider the following intermittent 
reduced map: 



R(x) = < 



x-3/4+c 

2c-l/2 
z-3/4+e 
l-2c 

x - 4 z ~ 1 (-a;) 2 
x + 4 2 " 1 - * 

x — c 
2c-l/2 
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5,1/2]. 



(17) 



The invariant density of the reduced map is given by 
dfi = h(x)\x\ 1 ~ z dx, which means an infinite invariant 
measure for z > 2 [22[ . Since the observation function of 

1/2 

a drift, f m (x), is an L 1 {^jl) function and J Q f m {x)d\i ^ 
(c / 0.375), the distributional limit theorem, Eq. (fTJ, 
holds. In particular, the distribution of the normalized 
time average of f m {x) converges to a Mittag-Leffler dis- 
tribution: 



n-l 



n .. — * 



M a , 



(18) 



where the return sequence is given by 



logn 



a n oc 



Thus, 



(z = 2) 



(z>2). 



(19) 



AA 



g(x)dn 



g c (x)dn = AA, 7 c Ai (J c ), (14) 



Moreover, the observation function of the Lyapunov ex- 
ponent, g(x) = ln|T"(x) = ln|i?'(x)|, is also the L+(/j.) 
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FIG. 2. (Color online) Einstein relation and constant \v 
(z — 2.5). Circles show the results of numerical simulations. 
Green lines show the theoretical results obtained by Eq. (J6]) 
and 1)160 . Constants (D) and (A) are obtained by numerical 
simulations. K is a fitting parameter (K = 1.12). 
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FIG. 3. (Color online) Relation between the difference of the 
generalized Lyapunov exponent and the ensemble-averaged 
velocity (z — 2.5). Circles are the results of numerical simu- 
lations and the line is the theoretical curve 1)150 . where Knax 
and ^{J c ) are obtained by numerical simulations. In particu- 
lar, we calculate /i(J c ) by lim„^. oo (^^ =0 ljc(i? fe a;)/a n ), and 
we set a n = n a . 



function. Thus, the normalized Lyapunov exponent is 
intrinsically random and its distribution converges to a 
Mittag-Lcfflcr distribution. 

Using numerical simulations, we confirmed the Ein- 
stein relation and the constant xv-, which are in good 
agreement with theory (Fig. 2). Moreover, the relation 
between the difference of the generalized Lyapunov ex- 
ponent and the ensemble-averaged velocity is also valid 
except for a large V (Fig. 3). The generalized Lyapunov 
exponent is maximized at V = and decreased according 
to the increase in |V|. 

Conclusion. — We derive dynamical systems corre- 
sponding to CTRWs. In a biased model, we obtain the 
Einstein relation for the time-averaged velocity and the 
time-averaged diffusion coefficient. Moreover, the differ- 
ence of the generalized Lyapunov exponent between unbi- 



ased and biased dynamical systems is represented by the 
ensemble-averaged velocity. Using Hopf's ergodic theo- 
rem, we find that the ratio between the time-averaged 
velocity and the Lyapunov exponent converges to a uni- 
versal constant. The universal constant is proportional 
to bias and the proportional constant is given by the dif- 
fusion coefficient and the Lyapunov exponent without a 
bias. In general, the ensemble-averaged velocity is rep- 
resented by the probability p: V = (2p — l)V max , and 
p = /i(J+)//x( J c ), where = [c, 1/2]. The relation be- 
tween AA and V will be universal if the map on cell, 
[—1/2 + L, L+ 1/2], is continuous because AA is related 
to fi(J+)/[i(J c ). Moreover, when the derivative of the 
map on the interval representing bias is the same as that 
of the unbiased map, the constant K is almost unity. 
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